Rotational dynamics and friction in double-walled carbon nanotubes 
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We report a study of the rotational dynamics in double-walled nanotubes using molecular dy- 
namics simulations and a simple analytical model reproducing very well the observations. We show 
that the dynamic friction is linear in the angular velocity for a wide range of values. The molecular 
dynamics simulations show that for large enough systems the relaxation time takes a constant value 
depending only on the interlayer spacing and temperature. Moreover, the friction force increases 
linearly with contact area, and the relaxation time decreases with the temperature with a power 
law of exponent —1.53 ± 0.04. 

PACS numbers: 68.35.Af; 85.35.Kt 
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Since the pioneering work of Iijima on nanotubes 0, 
research on their electronic and mechanic properties has 
been rapidly growing. Very high electric and thermal 
conductivities were observed and the nanotubes were 
shown to be as stiff as diamond in their axial direction 
with Young modulus of the order of TPa. These observa- 
tions make the nanotubes promising candidates for nano- 
electrical and nano-mechanical devices. The experiments 
on the translational motion in multi-walled nanotubes of 
Cumings and Zettl and Yu et al. 3] motivated the 
study of these possible oscillators. Many numerical stud- 
ies were performed to characterize the translational mo- 
tion of double- walled nanotubes 0, 0, IS El • O n the 
other hand, experiments were also carried out to study 
the rotational motion of concentric nanotubes. Indeed, 
Fennimore et al. and Bourlon et al. built actua- 
tors based on multi- walled nanotubes. As for the transla- 
tional oscillatory motion, no wear or fatigue was observed 
in the rotation. Recently, Krai and Sadeghpour [l^ pro- 
posed to spin nanotubes with circularly polarized light, 
further motivating the use of nanotubes as possible axis 
of rotation at the nanoscale. However, few systematic 
analytical and numerical studies of rotation have been 
done, except the work of Zhang et al. [13 where nonequi- 
librium molecular dynamics simulations were carried out 
to calculate the energy dissipation rate during the ro- 
tational motion in the presence of external thermostat. 
Nevertheless, the friction between two carbon nanotubes 
which are rotating one with respect to the other has not 
yet been understood and characterized. 

The purpose of the present Letter is to study the ro- 
tational dynamics of a double-walled nanotube (DWNT) 
by molecular dynamics simulations with the conserva- 
tive Hamiltonian mechanics in order to avoid additional 
sources of dissipation and precisely characterize the ro- 
tational friction. Since friction is related to the fluctu- 
ations of the many vibrational degrees of freedom, we 
set up a Langevin-type stochastic model and validate it 
by molecular dynamics simulations of typical DWNTs. 
This shows that dynamic friction is proportional to the 



angular velocity up to high values, increases linearly with 
contact area, and as a power law of exponent 1.53 ± 0.04 
with temperature. 

We consider a DWNT with free ends in vacuum. The 
Hamiltonian of the system of two concentric nanotubes is 
the same as in our previous studies for the translational 
motion 0, The intratube interactions are modeled 
by the Tersoff-Brenner potential with the set of param- 
eters of Table III in Ref. 01 an d the intertube interac- 
tions by a 6-12 Lennard- Jones potential used for carbon 
material simulations 0|. The DWNT is long enough 
so that the inner and outer nanotubes essentially rotate 
around a common axis with the angular velocities ui\ 
and u>2, respectively. Their angular momenta are given 
by L a — I a uj a in terms of the moments of inertia I a 
(a = 1,2). The equations of motion of this model are 
given by 



dL a duj a 

— J - — Af 

dt * dt _ " 



(a = 1,2) 



(1) 



with the torques N\ = —N2 in order to conserve the total 
angular momentum Lq = Li + L%. We can introduce the 
angular velocity of the center of inertia Q = (Jx^i + 
/2W2) + I2) = Lq/(Ii + 12) and the relative angular 
velocity oj — lo\ — L02 so that the equations of motion 
become dw/dt = N±/I and dtt/dt = with the relative 
moment of inertia / = hh/ih+h)- The total rotational 
kinetic energy of the nanotubes is given by 



1 2 1 



1 9 

2^ 



(2) 



The corrugation against the rotational motion is very 
small compared to the thermal energy. Indeed, for 
the two DWNTs considered in this Letter, namely 
(4,4)@(9,9) and (7J3)@(9,9), and the 6-12 Lennard- Jones 
potential we use |lj], we find a corrugation against ro- 
tation of respectively 0.00992 meV/atom and 0.00548 
meV/atom. These values are of the same order as the 
ones of Zhang et al. who used the Kolmogorov and Crespi 
potential |l5j. Even for temperatures as low as 100 K 



2 



(8.6 meV), thermal energy is much larger than corruga- 
tion energy against rotation. This observation allows us 
to neglect the mean static torque due to the corrugation 
in the equation of motion of the system. However, in or- 
der to model the rotational dynamics of the DWNT one 
should take into account dynamic friction. If the angu- 
lar velocity is small enough, we expect that the friction 
torque is linear in the angular velocity in analogy with 
the case of the translational motion Indeed, the 

rotation of a nanotube inside another one is comparable 
to the the sliding friction of Xenon film on a silver surface 
considered by Daly and Krim 0] and Tomassone et al. 
|l7j | who showed that, in their system, friction is domi- 
nated by phonon excitations and that dynamic friction 
is proportional to the sliding velocity. Accordingly, we 
model our system by a Langevin-type equation without 
potential, 

I ^7 = -X"> + WfluctM ( 3 ) 
at 

where % is a damping coefficient and Afl uct (t) a fluctuat- 
ing torque of zero mean value and correlation function 
(AWtWAWt(i')) - 2 X k B TS(t - t') in terms of the 
temperature T and Boltzmann's constant /cb- This is 
justified because the correlation time is of the order of 
the inverse of the Debye frequency (about 50 fs) as for 
the translational motion 4] . Since we observe relaxation 
times of the order of the nanosecond, one can assume 
this Markovian form for the evolution equation. We no- 
tice that the variations of the moments of inertia are 
neglected in Eq.(j3J). 

A molecular dynamics simulation can be carried out 
by giving an initial angular momentum Lq — Iiujq to the 
inner tube (a = 1) and a zero one to the outer tube. 
After the relaxation time r = I /x, the system reaches a 
state of equilibrium in which both tubes rotate with the 
same angular frequency, (wi) cq = (o^cq- This numerical 
observation is in agreement with the consequence of Eq. 
© that the mean relative angular velocity vanishes at 
equilibrium. Using the conservation of total angular mo- 
mentum, the equilibrium values of the angular momenta 
are thus given by (L a ) cq = L Q I a /(h + J 2 ) for a = 1,2. 
At equilibrium, the fluctuations of the angular velocity 
have the variance (oj 2 ) eq = k^T jl. For times longer than 
the relaxation t I/Xi the relative angle 9 performs a 
diffusive random rotation 

(9 2 (t))c2D(t-r) (4) 

with the diffusion coefficient D = ksT/x- Moreover, the 
angular velocity autocorrelation function is given by 

(c(iV(t'))eq=^ e - A(t -*' ) (5) 

where we defined A = \jl = 1/t. If a relative angular 
velocity loq is initially given to the system, Eq. shows 



that it decays on average according to 

(u(t)) = lo e~ xt (6) 
and the total rotational kinetic energy J2J) as 

<T,(t)> = 7^7- {he-^ + h) + ^ (1 - e-^) (7) 

where To = (T r (0)) denotes its initial value. Hence, we 
conclude that, in the long-time limit t — > oo, the rota- 
tional kinetic energy reaches the value 

We now show the results of the molecular dynamics 
simulations. The integration is done with a velocity Ver- 
let scheme of time step 2 fs for temperatures up to 600 K 
and 1 fs for larger temperatures. Angular velocity is fixed 
by giving a tangential velocity corresponding to coq to the 
velocities of each atom of the inner tube. We compute 
the average relative angular velocity (u>) and rotational 
kinetic energy (T r ) from the angular momenta and mo- 
ments of inertia in the axial direction. The angle is then 
obtained by the time integral of the angular velocity. 

As a first test, we estimated the maximum angular 
velocity the nanotube could endure before its disloca- 
tion. One can easily understand that the velocity of ro- 
tation should be smaller than the radial phonon velocity 
to guarantee the stability of the system. Indeed, we have 
calculated the critical angular velocity at which the nan- 
otube collapses for increasing radii of the nanotube. We 
found out that the critical velocity v c = u c R is constant 
and approximately equal to 7.98±0.80 km/s, close to the 
vel ocity 8 km/s of radial phonons However, Zhang et 
al. |l2( observed the collapse of their nanotube for angu- 
lar velocities higher than 4.7 rad/ps, corresponding to a 
radial velocity of approximately 3.1 km/s. This smaller 
value is probably due to the fact that the angular veloc- 
ity was applied incrcmcntaly to constraint atoms at the 
end of the tube while we give as an initial condition the 
angular velocity to all the atoms of the nanotube. 

We apply our model to the armchair-armchair DWNT 
(4,4)@(9,9) with both nanotubes of length 2.11 nm and 
inertia moments l\ = 139.11 amu nm 2 and I2 — 1547.9 
amu nm 2 (7 = 127.64 amu nm 2 ) at room temperature 
(300 K). We depict in Fig. ^the time evolution of the 
angular frequency and the rotational kinetic energy from 
an initial relative angular velocity of u>o = 1.07 rad/ps. 
The relaxation time calculated from the angular velocity 
in Fig. ^a) is 473 ps while the relaxation time from the 
energy is 243 ps which is in very good agreement with 
the factor 2 appearing between the decay rates in Eqs. 
© and 0. The decay of the rotational kinetic energy 
is well fitted by Eq. JjJ with the initial kinetic energy 
To = 0.822 eV. On the other hand, the simulations show 
that the equilibrium kinetic energy takes the value 0.0816 
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FIG. 1: (Color online) Plots of the time evolution of (a) the 
average relative angular velocity {u) and (b) the average rota- 
tional kinetic energy (T r ) for the armchair-armchair DWNT 
described in the text. The dashed lines are the fits by (a) Eq. 
© and (b) Eq. £J. The averages are computed with about 
10 trajectories. 
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FIG. 2: (Color online) Plots of the time evolution of (a) the 
velocity autocorrelation function and (b) the mean square dis- 
placement of the angle for the armchair-armchair DWNT de- 
scribed in the text. The dashed lines are the fits by (a) Eq. 
10 and (b) Eq. The statistics are computed with about 
10 3 trajectories. 



eV, which is in excellent agreement with Eq. © of the 
model. In order to confirm the validity of the stochastic 
model described in this Letter, we use the same code, but 
with a vanishing initial angular velocity. The velocity au- 
tocorrelation function and the mean square displacement 
are calculated after a time of 1 ns longer than the relax- 
ation. The results are depicted in Fig. [2] The velocity 
autocorrelation function is well fitted by the exponential 
© with a relaxation time of 475.9 ps while the linear 
fit 011 the mean square displacement gives r = 471.3 
ps. These results are in excellent agreement with the pre- 
vious simulations. The small differences are due to the 



statistical error, and can be reduced with longer simula- 
tion times. These simulations successfully demonstrate 
that the dynamic friction against the rotational motion 
is proportional to angular velocity in DWNTs. 

We calculated the dependence of the damping coeffi- 
cient on contact area at room temperature (300 K) for 
increasing lengths of our DWNTs. In order to obtain 
this result, we measured the relaxation time r of the an- 
gular velocity and calculated the damping coefficient as 
X = I/t. The damping coefficient x is depicted as a 
function of the length of the tubes in Fig. [31 where we 
see that the linear fit is in excellent agreement with the 
simulations. The friction force does not vanish at zero 
contact area because of the attractive forces between the 
surfaces 0] . Having obtained the dependence on length, 
we investigated the dependence on radius. We calculated 
the relaxation time for increasing radii of zigzag-armchair 
DWNTs. Our numerical simulation gives no systematic 
dependence of the relaxation time on radius except the 
variations due to the different interlayer spacings con- 
sidered. Hence the dependence on radius of the damp- 
ing coefficient comes only from the inertia moment I. 
The inertia moment of a nanotube of length radius R 
and surface mass density a is very well approximated by 
/ ~ 2ir<r£R 3 . For the hexagonal lattice, the surface mass 
density is a = 4m/(3v / 3flcc) wnere m — 12 amu is the 
atomic mass of carbon and acc = 1-42 A the carbon- 

n 2 

carbon bond length. Hence, a = 4.55 amu/A . In order 
to recover the experimental conditions 0, llOj where the 
inner ring is fixed and the outer free, we suppose that 
I\ — > oo so that the relative inertia moment is I = 1%. 
Since the friction torque is N = xcj and the angular ve- 
locity u) = v/R, the friction force is given by 



a A 



t = = — ;r V = 

R R? t 



(9) 



with the contact area A = 2tt£R. For nanotubes long 
enough £ the relaxation time becomes independent 

of the length as shown in Fig. |3| and the shear stress 
can be written as F/A — uv/too. As an example, the 
shear stress undergone by a nanotube of radius 10 nm 
with a rotational frequency of 1 GHz is approximately 
0.047 MPa. This confirms the expectation that the fric- 
tion force increases linearly with the contact area A . 
Nevertheless, we notice that friction could be enhanced 
because of the presence of impurities and deformations 
due to the plate fixed to the outer tube in the experi- 
ments of Refs. 

We have also investigated the dependence of the damp- 
ing coefficient on the temperature by the method of the 
autocorrelation function. A smaller DWNT have been 
used to gain computation time, namely a (4,4)@(9,9) of 
length 9.8 A. Indeed, the relaxation time is very large 
for low temperature. Since the dynamic friction is due 
to phonon interaction in the DWNT, one expects it to 
increase with temperature. On the other hand, at zero 
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FIG. 3: (Color online) Dependence of the damping coefficient 
X on length. The circles and squares are the results of the 
simulations for respectively the armchair-armchair DWNT 
(4,4)0(9,9) and the zigzag-armchair DWNT (7,0)@(9,9). The 
linear fits have the same slope 0.044 amu nm/ps. The inset 
is the relaxation time versus length I, which is well fitted by 
r = Tool/ '(I + la). Too = 1-22 ns and Iq = 3.6 nm for the 
armchair-armchair DWNT. Too = 1-30 ns and Iq — 1.0 nm for 
the zigzag-armchair DWNT. 
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FIG. 4: (Color online) Dependence of the damping coefficient 
X on temperature. The circles are the simulation results and 
the line is a power law fit of exponent 1.53 ± 0.04. The inset 
is the relaxation time versus temperature. 



temperature, the rotational motion should become bal- 
listic and the diffusion coefficient D = k^T/x should 
diverge. Indeed, our results confirm these expectations 
since we find that the damping coefficient depends on 
the temperature as the power law x ~ T v with the ex- 
ponent v = 1.53 ± 0.04, as seen in Fig. 0] Accordingly, 



the relaxation time for long enough DWNTs behaves as 
Tqc, = t*(T*/T) v where the time r* depends mainly on the 
interlayer spacing of the DWNT and T* is some reference 
temperature. We notice that dynamic friction may be 
non- vanishing at very low temperature where it is domi- 
nated by electronic interactions and impurities. 

In conclusion, we proposed a stochastic model describ- 
ing very well the rotational dynamics of DWNTs. We 
showed that dynamical friction is proportional to the an- 
gular velocity if this later is not too high. Moreover, for 
long enough nanotubes, the friction force is proportional 
to the contact area and can be characterized in terms of 
a relaxation time which depends on the temperature 
and the interlayer spacing. For higher velocities, non- 
Markovian and nonlinear effects can arise resulting into 
an increased dissipation. 
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